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The Membership and Distance of the Open Cluster ColUnder 419 
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ABSTRACT 

The young open cluster Collinder 419 surrounds the massive O star, 
HD 193322, that is itself a remarkable multiple star system containing at least 
four components. Here we present a discussion of the cluster distance based upon 
new spectral classifications of the brighter members, UBV photometry, and an 
analysis of astrometric and photometric data from the UCAC3 and 2MASS cat- 
alogs. We determine an average cluster reddening of E{B — V) = 0.37 ± 0.05 
mag and a cluster distance of 741 ± 36 pc. The cluster probably contains some 
very young stars that may include a reddened M3 III star, IRAS 20161+4035. 

Subject headings: stars: individual (HD 193322, IRAS 20161+4035) — stars: 
early-type — open clusters and associations: individual (Collinder 419) 
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1. Introduction 

Collinder 419 is an open cluster in the constellation Cygnus. It was first designated 
as Barnard 794 (Barnard 1927), but the accepted name comes from the work of Collinder 
(1931) who classified the cluster as a /x Normae cluster. This is an intermediate type of 
grouping between multiple star systems and open clusters. The prototype cluster NGC 6169 
contains the bright star // Normae = HD 149038 and a host of surrounding fainter stars. In 
the case of Collinder 419, the bright central star is HD 193322 = HR 7767 = HIP 100069 = 
WDS 20181+4044. HD 193322 is a multiple star system with a close visual companion STF 
2666 AB (separation of 2'.'9 and probably orbitally bound) and two more distant companions, 
STF 2666 AC (34'.'0) and TAR 5 AD (49'.'6), which, if bound have very long orbital periods. 
McAlister et al. (1987) discovered that the A component is a speckle binary, CHR 96 Aa,Ab. 
An approximately 31 year orbital period for this pair was computed by Hartkopf et al. 
(1993). In addition, one of the components of the speckle binary is a spectroscopic binary 
with a 311.03±0.25 day period (FuUerton 1990; McKibben et al. 1998). FuUerton (1990) also 
obtained a spectrum of the B component that suggests that it may also be a spectroscopic 
binary. Thus, the central AB pair may contain as many as five individual stars. 

The primary of the system, Aal, has a spectral type of 09 V:((n)) (Walborn 1972), with 
the suffix indicating shghtly broadened fines. McKibben et al. (1998) tentatively identify 
the close spectroscopic companion Aa2 and the B component all as B-type stars. The 
magnitude difference of the speckle pair Aa,Ab is AV = 1.2 ±0.5 (Mason et al. 2009) which 
suggests that Ab is also an early B-type star. Both the Aa2 and Ab companions offer us the 
means to determine mass for the O star component Aal through measurement of the orbital 
motion. By themselves, neither an orbit of the Aa,Ab pair nor the spectroscopic Aal,Aa2 
pair will yield a mass, but, in conjunction with a distance determination, a mass can be 
determined. The Hipparcos Catalogue (Ferryman et al. 1997) has become the standard 
source for parallaxes, but it may be unreliable for some double stars (Shatskii & Tokovinin 
1998) and also for distant O stars (Schroder et al. 2004). Until a new space-based parallax 
engine is available sometime in the next decade we are left with classical techniques as the 
best way to determine 0-star distances. 

We are currently pursuing efforts to determine the orbits of the central triple through 
speckle and optical long basefine interferometry (Turner et al. 2008). Here we present a 
new assessment of the distance to the cluster that we will use in the forthcoming orbit 
and mass determination. In §2 we present spectral classifications for the brighter members, 
and we show fits of their spectral energy distributions that provide estimates of reddening 
and angular size. We then describe in §3 an astrometric and photometric study of cluster 
membership based upon data in the UCAC3 catalog (Zacharias et al. 2010) and the 2MASS 
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catalog (Cutri et al. 2003). We also present new UBV photometry of the fainter cluster 
members in 54. Our results are summarized in §5. 



2. Spectroscopy and Reddening of the Brighter Stars 

We obtained spectra for five of the brighter stars in the cluster in order to determine 
spectral classifications and to estimate the intrinsic colors of the stars. One of us (MVM) 
obtained spectra of the four blue stars of the central multiple system using the Kitt Peak 
National Observatory (KPNO) 2.1 m telescope and Cassegrain focus, GoldCam CCD Spec- 
trograph. These spectra were made during 2005 October and November with the No. 47 
grating (831 grooves mm~^) in second order to obtain a resolving power X/6X ~ 3000. The 
observed wavelength range was 4050 - 4950 A, a region that includes numerous H Balmer, 
He I, and He II lines in O- and B-type spectra. The data were reduced using standard rou- 
tines in IRAF@. The flux rectified spectra were compared directly to spectra from the atlas of 
Valdes et al. (2004) in order to estimate the spectral classification based upon the strengths 
of the H Balmer, Ca II, He I, and Mg II lines. Our results are given in Table 1 for the B, 
C, and D components of HD 193322. See Figures [D - [3] for the comparision between the 
target stars and the template spectra. Our spectra for the combined flux of Aal, Aa2, and 
Ab agree with the 09 V:((n)) classification from Walborn (1972), which is listed in Table 1. 
These classifications all agree within a subtype with earlier results (Burnichon 1975; Hoffleit 
1982; Abt & Cardona 1983). 

The very red star IRAS 20161+4035 is the brightest object in the f^^-band in this 
vicinity, and we were curious to determine if indeed it is related to the cluster. One of us 
(EDG) obtained three spectra of the object with the KPNO 0.9 m Coude Feed Telescope in 
2004 October. These spectra cover the range from 6460 to 7140 A with a resolving power 
X/SX ~ 9700. The average spectrum is compared to several other cool star spectra from 
Valdes et al. (2004) (also made with the KPNO Coude Feed telescope) in Figure IH We 
made a digital comparison with all the M2 - M5 spectra in the Valdes et al. spectral atlas, 
and the best fit (based mainly on the appearance of the strong TiO bands in this spectral 
region) was the spectrum of the M3 II star HD 40239 (although there were a number of other 
M3 HI spectra which made an almost equally good fit). Figure H] shows the good overall 
match made with the M2 II and M3 HI giants compared to the dwarf and supergiant spectra. 



^IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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We also measured the Ha equivalent width in these and similar stars to help establish the 
luminosity class (Eaton 1995). We found that among stars of similar spectral type that 
the mean equivalent width of Ha rang ed from 0.76 ± 0.04 A for 2 main sequence stars, to 
1.59 ± 0.08 A for 4 giant stars, and up to 1.89 ± 0.04 A for 3 supergiant stars. The measured 
Ha equivalent width of IRAS 20161+4035, 1.38 A, clearly places the star among those in the 
giant luminosity class, and we adopted this assignment in Table 1. Note that the spectrum 
of IRAS 20161+4035 shows the strongest Li I A6707 feature among the spectra illustrated 
in Figure m 

We can use these spectral classifications to estimate the intrinsic colors of the targets, 
compare these with observed multi-color observations, and find the reddening E{B — V) in 
each case. However, the available near-IR 2MASS photometry for HD 193322 corresponds 
to the total fiux for the entire Aal, Aa2, Ab, and B complex. Thus, we need to make some 
assumptions about the relative fiux contributions of each component in each filter band to 
derive magnitudes and fiuxes for the individual stars. The magnitude differences of the AB 
pair were determined in the Johnson BVRI bands through adaptive optics imaging by ten 
Brummelaar et al. (2000). Unfortunately, much less is known about the A sub- components. 
Mason et al. (2009) found a magnitude difference of /W = 1.2 ± 0.5 mag for the Aa,Ab 
speckle pair, and McKibben et al. (1998) used a statistical method to estimate the mass ratio 
and magnitude difference AV = 1.3 mag for the Aal,Aa2 pair in the spectroscopic binary. 
These magnitude differences (plus AV for the AB pair from ten Brummelaar et al. 2000) 
lead to l^-band flux ratios relative to star Aal that are given in Table 2. To derive the flux 
ratios in other bands, we estimated their colors according to their spectral classiflcation and 
the color calibration from Wegner (1994). We assumed that the Aa2 and Ab stars are main 
sequence objects and then estimated their spectral types (B1.5 V and Bl V, respectively) 
based on their magnitude differences relative to Aal (using the absolute magnitude versus 
spectral type relation for main sequence stars from Lesh 1979). Table 2 lists the adopted 
flux ratio estimates for the sub-components (which agree with those from ten Brummelaar 
et al. for the Johnson R, I bands, but not for their i?-band result for which component B 
is unusually faint in B relative to the other bands). We used these flux ratios to determine 
the magnitudes of Aal and B, the stars with reliable spectral classiflcations, in the following 
analysis. Note that the current flux ratio uncertainties probably introduce an error of ±0.2 
mag into the error budget for the magnitudes of these two stars. 

We estimated the reddening and angular sizes of the flve bright stars listed in Table 1 
by comparing their observed fluxes with reddened model spectral energy distributions. The 
optical magnitudes were taken from the work of ten Brummelaar et al. (2000; for components 
Aal and B), Burnichon & Garnier (1976; for components C and D) and Droege et al. (2006; 
plus UBV magnitudes given in §4 for IRAS 20161+4035), and these were combined with 
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2MASS magnitudes for the near-IR (Cutri et al. 2003). The magnitudes were transformed 
to fluxes using cahbrations from Cohna et al. (1996; Johnson U,B,V), Bessell et al. (1998; 
Cousins Jc), and Cohen et al. (2003; 2MASS J, H, Ks). We used Kurucz model atmosphereqlj 
(assuming \ogg = 4.0, solar abundances, and a microturbulent velocity of 2 km s~^) for 
temperatures from the spectral calibration of Bohm-Vitense (1981; see Table 1) to create 
low resolution versions of the unreddened spectral energy distributions (SEDs). However, 
for the case of IRAS 20161+4035 (M3 III), we adopted a flux distribution from the MARCS 
cod (Gustafsson et al. 2008) for a model with T^s = 3500 K, logs' = 1-0, solar abundances, 
and a microturbulent velocity of 5 km s~^. We then fit the observed fluxes with model SEDs 
reddened according to the formulation from Fitzpatrick (1999). Because the stars have low 
reddening and limited short-wavelength coverage, it is difficult to obtain an independent 
estimate of the ratio of total-to-selective extinction R, so we adopted the default value 
i? = 3.1 throughout. The derived reddening E{B — V) and limb-darkened angular diameter 
6ld are given in Table 1. The blue stars yield a consistent reddening of E{B—V) = 0.29±0.01 
mag, which is significantly less than that found for IRAS 20161+4035, E{B-V) = 0.74+0.04 
mag. The SED for this red star is shown in Figure [5], and we see that there appears to be 
an IR flux-excess of 48 + 10% near 12/im, based upon data from the MSX (Egan & Price 
1996) and IRA^^ missions. We argue below (§5) that the higher reddening and IR-excess 
may be related to the youth of this star. 



3. Astrometric and Photometric Data from UCAC3 and 2MASS 

With the recent release of the third U. S. Naval Observatory CCD Astrograph CataloJ^ 
(UCAC3; Zacharias et al. 2010), we now have access to a large collection of astrometric and 
photometric data that can be used to search for additional cluster members. We began by 
selecting all the stars in UCAC3 with complete sets of proper motion and photometric data 
that are found within a 0?25 radius of HD 193322, the nominal center of Cr 419. We show 
in Figure [6] the average surface density of star counts within annuli centered on the position 
of HD 193322. The larger surface density associated with the cluster extends outwards to at 
least 0?10, and we set the outer boundary on cluster membership at 0?16 to avoid including 
too many field stars at large distance. Figure [7] shows a 0?25 field centered on HD 193322 



^"^httpi/Zkurucz. harvard.edu/ 
^^http://niarcs. astro. uu.se/ 

i^Infrared Astronomical Satellite Catalogs, 1988; The Point Source Catalog, version 2.0, NASA RP-1190 
^'^http://www. usno.navy.mil/USNO/astrometry/optical-IR-prod/ucac 
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and shows the low stellar density of the cluster. The image is from t he near-IR plates of 



Second Palomar Observatory Sky Survey (P0SS2-IR) (IReid et al.lll99ll ). Figure |8] shows the 
distribution of proper motions of stars found within this radius. The circle in the diagram 
is centered on fia cos 6 = —2.2 mas y~^ and fis = —9.1 mas y~^, the average values for the 
bright cluster stars (HD 193322A, C, D, and HD 228810). There appears to be a bimodal 
distribution in proper motion in this direction, so we selected only those stars within the 
circle, i.e., those with a combined proper motion difference from the cluster mean of less than 
11 mas y~^ (the RMS of the UCAC3 proper motion errors and the expected range within the 
cluster itself). The constraints of position and proper motion reduced the potential sample 
to a group of 96 stars. 

The histogram in Figure E] also demonstrates that there will be a significant number 
of field stars in this sample. To estimate this number, we applied the same proper motion 
constraint to the UCAC3 sample within a separation range of 0?16 to 0?24, and this indicates 
that the surface density of field stars in this direction with similar proper motions to the 
bright cluster stars is 1000 ± 100 stars per deg^. Thus, we might expect to find some 83 ± 8 
field stars in the 0?16 radius circle centered on HD 193322. This exercise suggests that the 
actual number of cluster stars is small (perhaps 96 total - 83 field = 13 cluster stars) and 
that the sample will be dominated by other field stars. 

We next considered the photometric data in UCAC3, which consists of 2MASS photom- 
etry and the UCAC3 red-band photometry that we will refer to as i?„ (to distinguish it from 
the Johnson and Cousins counterparts). This filter has a bandpass from 579 to 642 nm that 
falls between the Johnson V and Cousins Rc bands. The UCAC3 catalog reports magnitudes 
from centroid fit models and aperture photometry, and we used the latter in our analysis. 
The UCAC3 listing for HD 193322 corresponds to the total flux from components Aal, Aa2, 
Ab, and B, and we used the flux ratios given in Table 2 to determine the corresponding mag- 
nitudes for stars Aal and B only. Our goal is to compare the magnitudes of the candidate 
member stars with model predictions for given assumptions about the reddening E{B — V) 
and distance modulus DM, in order to find values of these parameters that lead to the best 
fit of the photometry. 

We adopted model isochrones for the cluster from the evolutionary models of Bertelli 
et al. (1994) and Marigo et al. (2OO8)0. We selected an isochrone for solar metallicity and 
a cluster age of 7 Myr, the value proposed by Kharchenko et al. (2005) and one that is 
consistent for a time of flight from Cr 419 of the runaway star 68 Cyg (Schilbach & Roser 
2008). The isochrones were produced for a number of photometric bands including Johnson 



14 



http://stev.oapd.inaf.it/cmd 
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V, Cousins and 2MASS J, H, K^. 

We need to transform the model optical magnitudes into the Ru system in order to 
compare them with the UCAC3 values. We used a sample of 402 photometric standard stars 
from Landolt (1992) that appear in UCAC3 to arrive at a transformation equation from 
Johnson V and Cousins Rc to UCAC3 Ru, 

Ru = Rc + 0.292 + 0.372(F - + 0.054(i?c - 13) (1) 

which successfully fits the observed UCAC3 Ru magnitudes over the range i?„ = 8 — 17 
with a standard deviation of 0.116 mag. The color term in the equation shows how the Ru 
magnitude changes with stellar color because of the difference in the effective wavelengths 
of the Ru and Rc bands. The final term on the right hand side reflects a small systematic 
difference between the Landolt and UCAC3 systems that is approximately equal to the 
scatter in the relation at the bright and faint ends of the sample. 

We can apply a few other constraints to limit the number of field stars that may exist 
within the astrometric group of cluster candidates. Figure M shows the near-IR color-color 
diagram for the sample together with unreddened {solid line) and reddened {dashed line) 
versions of the cluster isochrone. The numbers beside the unreddened version indicate the 
stellar mass corresponding to the color pair, and massive stars have colors in the lower part of 
the diagram. The reddening vector direction is indicated by the dotted line. Comeron et al. 
(2002) and others have argued that the reddening-free parameter Q = {J— H) — 1.9 A* {H— Kg) 
(Fitzpatrick 1999; Indebetouw et al. 2005) provides an excellent means to identify massive, 
reddened stars that typically have Q ~ 0.0. We decided to apply this criterion by restricting 
the sample to those massive stars with Q < 0.116, the model value for a star of mass 
M = 1.4M0. Imposing this limit will mean that we lose fainter, red cluster members, 
but, on the other hand, we will also remove cool background giants and cooler, low mass, 
foreground stars that do not belong to the cluster (for example, those stars found near the 
unreddened isochrone between 0.8 and I.4M0). 

We also applied a faint limit constraint, Ru < 15, to avoid stars with large magnitude 
errors. Since our initial sample was drawn from UCAC3 stars with complete photometric 
coverage, the result is that this sample tends towards stars with equal magnitude in all filters 
at the faint end. This group probably includes foreground A-type stars, so applying the faint 
limit will help us remove such stars from consideration. 

We show two versions of the color - magnitude diagram in Figures [TU] & [TTJ Both 
of these show evidence of a main sequence group at the blue end, but both also show a 
large population of redder stars that lie above the expected main sequence. In order to 
remove such targets, we restricted the sample to those with i?„ — J < 1.35. This constraint. 
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together with the Q hmit described above, will result in a sample without any low mass, 
main sequence members. 

Our final constraint is based on how closely the magnitudes match the model isochrone 
values. For any given choice of E{B — V) and DM, we formed an average statistic between 
all the observed and model magnitudes (reddened according to the relations given by Fitz- 
patrick 1999) for each mass point along the model isochrone. The best fit (lowest x^) model 
point was taken as an appropriate mass for the specific star. Those stars with magnitudes 
and colors very different from the main sequence values will have a large minimum x^, and 
we restricted our final sample of main sequence cluster members to those with < 35, 
a value high enough to include those cases where the observational photometric errors are 
underestimated. The final sample consists of 16 stars that are listed in Table 3 and are 
indicated in Figures 191 - ITT] by plus signs. Note that we make the tacit assumption that the 
same reddening can be applied to all cluster members, and this can only be true in some 
average sense. There may be a reddening gradient across this field (Schlegel et al. 1998) and 
some of the youngest objects may suffer from circumstellar dust reddening. Nevertheless, 
until spectroscopy is available for many stars and their intrinsic colors found, we must accept 
this working hypothesis. 

We performed a grid search over a range of E{B — V) and DM to find those values that 
minimized the average value for the selected cluster stars. The best fit was made with 
E{B — V) = 0.37 ± 0.05 mag and DM = 9.35 ± 0.03 mag, yielding an average photometric 
statistic of x^ = 4.4. The isochrones for these parameters are plotted in Figures [9] - [H], and 
the corresponding main sequence masses from the isochrone are listed in the final column of 
Table 3. Our derived reddening estimate is mainly consistent with earlier estimates for the 
cluster and star {E{B -V) = 0.34, Kharchenko et al. 2005; E{B -V) = 0.345, Burnichon 
1975; E{B -V)= 0.38, Schr5der et al. 2004; E{B -V) = 0.41, Cardelh et al. 1989). The 
reddening is slightly larger than that found by SED fitting for the blue stars (Table 1), but 
we suspect that the differences arise in the variable reddening across the cluster. We also 
repeated the analysis using an isochrone for a younger cluster (age 3 Myr), but the results 
were unchanged because the predicted colors and magnitudes were only significantly different 
for the most massive star HD 193322Aal (the predicted colors are marginally closer to the 
observed ones for the 7 Myr isochrone). 

There are many stars that are plotted above the main sequence in Figures [10] & [11] 
Some of these stars may belong to a massive and more distant population. For example, in 
the infrared color-color diagram (Fig. [9]), there appears to be a reddened group of OB stars 
near H — Kg = 0.2 and J — H = 0.4. An independent analysis of the 2MASS magnitudes of 
stars in the vicinity of Cr 419 by David Turner (priv. communication) suggests that these 
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correspond to objects with E{B — V) = 1.5 at a distance of 1.6 kpc that niay be associated 
with the nearby cluster Berkeley 87 (ITurner fc Forbes! Il982l : iMassey et al.ll200ll ). We show 
how the 7 Myr isochrone would appear for such reddened and distant stars as a dotted line 
in Figures nH] and [TT| and many of the stars appearing above the Cr 419 main sequence are 
close to the predicted colors and magnitudes of a distant population. Cr 419 is probably 
young enough that it may host pre-main sequence stars that would also appear above the 
main sequence. Dashed lines in Figure TOl and [TT] indicat e the pre-main sequence isochrone 
for an age of 7 Myr from the work of ISiess et al.l (120001 ) transformed to the distance and 



reddening of Cr 419, and indeed we find many examples of stars close to the predicted track. 
We suspect that redder stars appearing in these color-magnitude diagrams probably include 
some cluster pre-main sequence stars and large numbers of distant massive stars and other 
line-of-sight field stars. 



4. UBV Photometry 

We also obtained new UBV photometry of the central region of the cluster to explore 
the optical colors of the fainter stars. One of us (JRP) collected CCD images of Collinder 419 
in 2006 October with the Emory University Observatory, DFM 0.6 m Cassegrain telescopj^. 
The detector was an Apogee Ap47 1024 x 1024 pixel CCD camera. In order to concentrate 
the signal for fainter sources, the camera was set to a 2 x 2 binning mode with a net pixel 
scale of I'.'l pixel"^ and a 9'.5 x 9^5 field of view. The images of the brightest stars were 
saturated in individual frames in order to obtain a stronger signal on the faintest sources. 
Thirty frames were collected through each Johnson UBV filter. In addition, standard stars 
from charts 140, 141, and 148 of Landolt (1992) were observed over a large range in air mass 
to transform the instrumental magnitudes to the standard system. 

All the frames were debiased, dark subtracted, and flat fielded using standard routines 
in IRAF. The images were combined using IMCENTROID and IMCOMBINE to build a 
master image in each filter. Photometric measurements of the standard stars were made 
with aperture photometry routines in IRAF. The transformation coefficients were computed 
using the IDL procedure transfl.pro, written by Marc W. Buie. This procedure uses the 
general formula 

mo = mi- k'X - k"CX + eC + Z 

where mo is the standard magnitude, rrij is the instrumental magnitude, k' is the first order 
extinction coefficient, X is the airmass, k" is the second order extinction coefficient, C is a 



^^http://www. physics. emory.edu/astronomy/observatory.html 
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color index, e is a color coefficient, and Z is the magnitude zero point. The instrumental 
magnitudes of the stars in the CoUinder 419 frames were measured using PSF fitting routines 
in IRAF, and the relation between the PSF and aperture measurements was determined for 
ten well isolated stars in the field. The transformation equations were then used to determine 
Johnson magnitudes for the cluster stars. The errors in the transformation coefficients were 
added in quadrature with the instrumental magnitude errors to obtain the net observational 
errors (approximately 0.04, 0.03, and 0.02 mag in UBV, respectively). The results are listed 
in Table 4. 

We show in Figure [12] the color - magnitude diagram {B — V, V) for 43 of the 79 stars 
measured that meet the proper motion limits for cluster membership described above (§3). 
Most of the proposed cluster members (§3; Table 3) were too bright and/or outside of the field 
of view for our U BV photometric measurements, and we have U BV photometry for only 5 of 
the 16 proposed members (indicated by plus signs in Fig. 8 and by "Yes" in the final column 
of Table 4). Figure [T2l also shows the main sequence relation from the models of Bertelli et al. 
(1994) and Marigo et al. (2008) for our derived estimates of E{B -V) = 0.37 mag and DM 
= 9.35 mag, and this relation appears to be consistent with our optical photometry. There 
are a number of stars that have magnitudes indicative of lower mass, main sequence stars, 
but there are many objects with a position above the main sequence with a distribution 
similar to that seen in the red and near-infrared color - magnitude diagrams (Figs. [10] and 
[TT]) . Figure [13] shows the optical color - color diagram {B ~ V,U — B) for the subset with 
complete UBV measurements and within the proper motion limits for membership. We also 
show both unreddened and reddened versions of the model colors. Again, there is evidence 
for some lower main sequence members (towards the middle and lower sections of the dashed 
line), but there is also a group of stars with colors near {B — V,U — B) ^ (1.0, 0.0) that are 
characteristic of B-stars with a reddening E{B — V) ^ 1.0 mag. Most of these stars also 
have positions in the near-IR color - color diagram associated with reddened, intermediate 
mass stars (Fig. U{H -Ks,J-H)^ (0.15, 0.30)). 



5. Discussion 

The distance we obtained from the photometric analysis of d = 741 ± 11 pc (formal 
error) does not account for possible systematic errors in the model isochrone magnitudes. 
By comparing the {Rc — J, J) color- magnitude diagrams for a 7 Myr cluster from Marigo et al. 
(2008) with that from Lejeune & Schaerer (2001), we estimate that the model uncertainties 
amount to A J ^ 0.10 mag over the main sequence range of interest, so we add this error in 
quadrature with the fitting errors to obtain d = 741 ±36 pc. We can check this with distances 
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derived from the angular sizes given in Table 1. This distance estimate depends on the ratio 
of the physical diameter to the angular diameter, and for this comparison, we estimated the 
stellar radii for HD 193322B, C, and D from the model isochrone values for main sequence 
stars with our adopted effective temperatures (Table 2). We could not use this relation for 
HD 193322Aal because the model isochrone did not extend to stars quite this hot, so we 
instead adopted a radius of = 7.53Rq from the spectral calibration for main sequence 
stars from Martins et al. (2005). We estimated a radius of i? = 75Rq for the MS III star 
IRAS 20161+4035 from the cahbration for giants by van Belle et al. (1999). The resulting 
distances are hsted in the final column of Table 1, and the average, d — 766 ± 51 pc, agrees 
with the photometric result. Our results confirm the only previous distance estimate for the 
cluster from Kharchenko et al. (2005) of 740 pc (no error quoted). The cluster distance also 
agrees within errors with the Hipparcos distance for HD 193322 of 600lj^iQ pc (van Leeuwen 
2007). 

We suspect that some of the red stars that meet the spatial and proper motion con- 
straints of cluster membership but lie above the main sequence are very young objects. A 
similar population of bright, red objects has been reported for other young star clusters 
containing O stars (Kumar et al. 2004; Comeron & Pasquali 2005; Carlson et al. 2007; 
Negueuela et al. 2007; Comeron et al. 2008). These stars may include IR-excess objects with 
remnant disks like the Herbig Ae/Be stars. Although it is possible that the red object, IRAS 
20161+4035, is an old star that is coincidentally passing through the vicinity of Cr 419, we 
suspect that it may also be a very young member of the cluster, despite its evolved appear- 
ance. The star's large reddening and IR-excess both suggest that it is immersed in natal 
dust, and its strong Li I A6707 and high luminosity may be indicative of youth. Kumar et 
al. (2004) identified three other cases of luminous, late-type giants that are found in young 
clusters. These cool giants may represent another kind of spectroscopic manifestation of 
young stellar objects. 
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Fig. 1. — Our estimate f o r the sp ectral type of HP 19332 2 B is B 1.5 V based on comparisons 



to stars in 



Valdes et al.l ( 120041 ). IWalborn fc Fitzpatrickl ( Il990l ) advocate using the Si lines 



to classify early B stars, and those lines were given the strongest weight. The spectra were 
normalized and then offset to improve clarity. 
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Fig. 2.— We estimate the spectral type of HD 193322 C as B8 V or possibly B9 V. The 
spectra were normalized and then offset to improve clarity. 
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Fig. 3. — We estimate the spectral type of HD 193322 D as B9 V due to the strength of the 
Ca II K hnes. The H Balmer, He, and metal lines are weak. The spectra were normalized 
and then offset to improve clarity. 
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Fig. 4. — A comparison of the red spectrum of IRAS 20161+4035 with similar spectra of 
cool stars from the atlas of Valdes et al. (2004). Each spectrum has been offset in relative 
intensity for clarity of presentation. 
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Fig. 5. — The spectral energy distribution of IRAS20161+4035 {plus signs) compared with 
a reddened model SED {solid line). 
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Fig. 6. — The distribution of the stellar number density per unit area in rings of increasing 
radius centered on HD 193322, the nominal center of Cr 419. 
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Fig. 7.— The P0SS2-IR image centered on HD 193322 with a field of view of a 0.25° x 0.25° 
showing the low stellar density of the cluster. North is up and east is to the left. 
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Fig. 8. — The distribution of proper motions in right ascension (/x^ cos S) and dechnation 
(lis) for stars within 0?16 of the cluster center. The position of the circle marks the cluster 
average proper motion and its circumference represents the boundary for possible cluster 
membership. 



-24- 




Fig. 9. — The distribution of near-IR colors for those stars meeting the positional and proper 
motion criteria for cluster membership. Plus signs mark those that are probable, main 
sequence members while x signs are used for the rest of the sample. The solid and dashed 
lines represent the values for the cluster isochrone for no reddening and E{B — V) = 0.37, 
respectively. The numbers beside the unreddened isochrone indicate the stellar masses at 
the terminal points and at the origin of the reddening trajectory {dotted line). Stars below 
the dotted line are candidate massive stars (M > IAMq). 
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Fig. 10. — A color- magnitude diagram based upon the UCAC3 i?„ and 2MASS J magni- 
tudes. Plus signs mark probable, main sequence members, while x symbols are used for 
the remainder of the sample that meets the positional and proper motion criteria for cluster 
membership. The solid hne shows the best fit isochrone. The dashed line indicates the 
pre-main sequence isochrone for an age of 7 Myr. The dotted line is the isochrone for distant 
OB stars that may be related to the nearby cluster Berkeley 87. 
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Fig. 11. — A color-magnitude diagram for 2MASS J and Kg magnitudes (in the same format 
as Fig. [ID]). The bright, red object {top, right) is IRAS 20161+4035. 



-27- 




B - V (mag) 



Fig. 12. — A color-magnitude diagram from optical photometry in the same format as 
Fig. [TUJ The solid line shows the isochrone for our derived E{B — V) and DM. 
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Fig. 13. — A color-color diagram from optical photometry in a format similar to Fig. |9l The 
sohd and dashed lines show the model unreddened and reddened colors, respectively. 
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Table 1. Reddening Estimates for i? = 3.1 



Star 


Spectral 


E{B - V) 


dLD 


d 


Name 


Classification 


(mag) 


{fias) 


(kpc) 


HD 193322Aal . . . 


09 V:{(n)) 


0.31 ±0.02 


85 ±2 


0.83 


HD 193322B 


B1.5 V 


0.28 ± 0.03 


53 ±2 


0.73 


HD 193322C 


B8 V 


0.29 ± 0.02 


23 ± 1 


0.79 


HD 193322D 


B9 V 


0.29 ± 0.03 


25 ± 1 


0.70 


IRAS 20161+4035 


M3 in 


0.74 ± 0.04 


893 ± 46 


0.78 



Table 2. Adopted Flux Ratios for the Components of HD 193322 



O Ldl J. Qff 




Fx/Fj^al 


Fx/F^al 


Fx/FAal 






]\.Tamr> (lrK\ 

i\ ciriio (^Krv J 


B 


V 


Rj 


h 


J 




H 




HD 193322Aal 33.2 


1.00 


1.00 


1.00 


1.00 


1.00 




1.00 


1.00 


HD 193322Aa2 23.0 


0.29 


0.30 


0.31 


0.33 


0.33 




0.35 


0.35 


HD 193322Ab. 25.0 


0.41 


0.43 


0.44 


0.46 


0.47 




0.48 


0.48 


HD 193322B . . 23.0 


0.21 


0.21 


0.21 


0.21 


0.23 




0.24 


0.25 




Table 3. 


Main Sequence 


Cluster Members 








UCAC3 


Other 


Ru 


Ru — J 


J-H 






Mass 




Number 


Name 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


{Mq) 




262-202014 




12.29 


1.06 


0.14 


0.08 


11.02 


2.2 




262-202033 




12.90 


1.20 


0.23 


0.12 


11.35 


1.7 




262-202070 


HD 228810 


10.01 


0.55 


0.10 


0.05 


9.31 


5.0 




262-202090 




12.44 


0.94 


0.08 


0.09 


11.34 


1.9 




262-202099 




13.66 


1.33 


0.26 


0.09 


11.98 


1.4 




262-202103 




12.91 


1.14 


0.21 


0.12 


11.44 


1.7 




262-202123 


HD 193322C 


11.14 


0.51 


0.02 


-0.02 


10.63 


2.8 




262-202132 


HD 193322D 


11.26 


0.58 


0.03 


0.08 


10.57 


2.8 




262-202133 


HD 193322Aal 


6.59 


0.07 


-0.07 


-0.04 


6.64 


14.9 




262-202133 


HD 193322B 


8.29 


0.17 


-0.03 


0.01 


8.14 


9.0 




262-202155 




12.27 


0.71 


0.07 


0.02 


11.46 


1.8 




262-202269 




12.25 


0.72 


0.06 


0.07 


11.40 


1.8 




262-202314 




12.83 


1.09 


0.23 


0.06 


11.45 


1.7 




262-202350 




12.37 


0.88 


0.06 


0.10 


11.34 


1.9 




262-202365 




11.63 


0.74 


0.07 


0.05 


10.77 


2.5 




262-202402 




11.67 


0.78 


0.11 


0.02 


10.76 


2.5 





Table 4. UBV Photometry 



UCAC3 


U 


B 


V 




Number 


(mag) 


(mag) 


(mag) 


Member 



262- 


■201944 


12.51±0.04 


12 


.68±0 


.02 


11, 


.58±0, 


.01 


262- 


■201945 


14.26±0.04 


14, 


.42±0, 


.02 


13, 


.59±0, 


.02 


262- 


■201948 




15, 


.92±0, 


.03 


13, 


.91±0, 


.02 


262- 


■201953 




16 


.52±0, 


.04 


15. 


.21±0. 


.01 


262- 


■201963 


16.25±0.05 


16, 


.09±0, 


.03 


14. 


.91±0. 


.02 


262- 


■201973 


16.22±0.06 


16, 


.13±0, 


.04 


15. 


.01±0. 


.02 


262- 


■201978 


16.60±0.06 


16 


.llitO 


.04 


15, 


.02±0, 


.02 


262- 


■201981 


15.30±0.04 


15, 


.27±0, 


.02 


14, 


.43±0, 


.01 


262- 


■201993 


16.04±0.05 


15 


.76±0 


.02 


14, 


.83±0, 


.02 


262- 


■201997 




16 


.66±0 


.03 


15, 


.53±0, 


.03 


262- 


■202009 




17, 


.15±0, 


.05 


15. 


.86±0. 


.01 


262- 


■202014 


12.90±0.04 


12 


.83±0 


.01 


12, 


.41±0, 


.03 


262- 


■202017 


15.74±0.04 


15 


.23±0 


.02 


14, 


.18±0, 


.01 


262- 


■202018 




14 


.99±0 


.04 


14, 


.67±0, 


.03 


262- 


■202025 




16, 


.47±0, 


.03 


14, 


.62±0, 


.01 


262- 


■202026 




16, 


.31±0, 


.03 


15, 


.22±0, 


.01 


262- 


■202027 




16, 


.49±0, 


.03 


15. 


.61±0. 


.01 


262- 


■202032 


14.72±0.04 


14, 


.64±0, 


.01 


13. 


.74±0. 


.01 


262- 


■202033 


13.65±0.04 


13 


.59±0 


.01 


13, 


.08±0, 


.02 


262- 


■202034 




16, 


.12±0, 


.02 


14, 


.04±0, 


.02 


262- 


■202039 


15.17±0.04 


14, 


.92±0, 


.01 


14, 


.08±0, 


.01 


262- 


■202042 




15, 


.95±0, 


.02 


15, 


.32±0, 


.02 


262- 


■202045 


14.40±0.04 


14, 


.20±0, 


.03 


13, 


.55±0, 


.03 


262- 


■202046 




16, 


.88±0, 


.04 


15. 


.73±0. 


.01 


262- 


■202047 




17 


.26±0, 


.05 


15. 


.96±0. 


.01 


262- 


■202050 




16 


.52±0 


.03 


15, 


.22±0, 


.01 


262- 


■202052 


15.33±0.04 


14, 


.99±0, 


.02 


14, 


.20±0, 


.03 


262- 


■202054 


16.21±0.05 


15, 


.74±0, 


.02 


14, 


.76±0, 


.02 


262- 


■202059 


14.90±0.04 


13, 


.90±0, 


.01 


11, 


.51±0, 


.04 


262- 


■202061 




17, 


.47it0, 


.05 


16. 


.64±0. 


.01 


262- 


■202066 


14.98±0.04 


14 


.73±0 


.02 


13. 


.88±0. 


.03 


262- 


•202067 


14.64±0.04 


14, 


.46±0, 


.01 


13, 


.75±0, 


.02 


262- 


■202070 


10.19±0.04 


10, 


.57±0, 


.01 


10, 


.33±0, 


.01 


262- 


■202088 


13.09±0.04 


13 


.03±0, 


.01 


12, 


.48±0, 


.02 


262- 


■202089 


16.76±0.07 


16, 


.24±0, 


.03 


15, 


.13±0, 


.02 


262- 


■202097 




16, 


.21±0, 


.03 


14. 


.60±0. 


.01 


262- 


■202098 


15.69±0.05 


15, 


.35±0, 


.02 


14. 


.53±0. 


.01 


262- 


■202100 




17 


.24±0 


.05 


15, 


.92±0, 


.02 


262- 


■202107 


16.02±0.05 


15, 


.67±0, 


.02 


14, 


.95±0, 


.02 


262- 


■202109 




16, 


.33±0, 


.03 


15, 


.31±0, 


.02 


262- 


■202117 


13.98±0.04 


13, 


.95±0, 


.01 


13, 


.00±0, 


.02 


262- 


■202118 




16, 


.08±0, 


.05 


15, 


.53±0, 


.01 


262- 


■202120 




15, 


.26±0, 


.02 


14. 


.11±0. 


.02 


262- 


■202128 


13.56±0.04 


13, 


.56±0, 


.01 


12. 


.62±0. 


.01 


262- 


■202135 


16.41±0.05 


16, 


.15±0, 


.02 


15. 


.12±0. 


.01 



Yes 



Yes 



Yes 
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Table 4 — Continued 



UCAC3 
Number 


U 
(mag) 




B 
(mag) 




V 
(mag) 




Member 


262-202141 


16 


66±0 


09 


15.72±0 


03 


14.48±0 


01 




262-202146 








16.35±0 


03 


15.13±0 


04 




262-202155 


12 


62±0 


04 


12.67±0 


01 


12.39±0 


03 


Yes 


262-202156 








16.79it0 


03 


15.72±0 


02 




262-202159 


16 


45±0 


05 


16.09±0 


02 


15.08±0 


01 




262-202161 








17.00±0 


03 


16.07±0 


02 




262-202166 


14 


60±0 


04 


14.44±0 


02 


13.88±0 


02 




262-202170 


11 


59±0 


04 


11.63±0 


01 


10.68±0 


03 




262-202175 








17.59±0 


08 


16.33±0 


05 




262-202176 


12 


57±0 


04 


12.61±0 


01 


11.64±0 


01 




262-202177 


1 

io 




04 


14.70±0 


02 


13.65±0 


04 




262-202178 








16.73±0 


04 


15.81±0 


01 




262-202185 








17.45±0 


06 


16.09±0 


01 




262-202187 








17.41±0 


06 


1d.30±0 


01 




262-202190 


14 


87±0 


04 


14.81±0 


01 


13.80±0 


01 




262-202194 


15 


97±0 


05 


15.68±0 


02 


14.71±0 


03 




262-202197 








16.89±0 


04 


15.79±0 


01 




262-202199 


14 


34±0 


04 


14.18±0 


01 


13.13±0 


01 




262-202207 


14 


56±0 


04 


14.39±0 


01 


13.68±0 


01 




262-202212 


15 


56±0 


04 


15.41±0 


02 


14.32±0 


02 




262-202218 


13 


14±0 


04 


13.21±0 


01 


12.33±0 


01 




262-202225 








16.96±0 


04 


15.83±0 


02 




262-202224 








16.75±0 


03 


15.66±0 


01 




262-202227 


16 


31±0 


05 


15.84±0 


02 


14.89±0 


06 




262-202236 


16 


47±0 


05 


16.18±0 


03 


15.41±0 


01 




262-202239 








17.00±0 


04 


15.81±0 


04 




262-202243 


16 


16±0 


05 


15.65±0 


02 


14.80±0 


01 




262-202247 


13 


70±0 


04 


13.72±0 


01 


12.85±0 


03 




262-202248 


13 


57±0 


04 


13.34±0 


02 


12.62±0 


07 




262-202255 








16.40±0 


03 


15.56±0 


01 




262-202256 








17.86±0 


10 


16.15±0 


03 




262-202260 








16.10±0 


03 


14.84±0 


01 




262-202265 


14 


56±0 


04 


14.08±0 


02 


13.27±0 


01 




262-202269 


12 


59±0 


04 


12.67±0 


02 


12.38±0 


02 


Yes 



